Pharmacophore-Based Virtual Screening, Quantum Mechanics Calculations, and Molecular Dynamics Simulation Approaches Identified Potential Natural Antiviral Drug Candidates against MERS-CoV S1-NTD.

Middle East respiratory syndrome coronavirus (MERS-CoV) is a highly infectious zoonotic virus first reported into the human population in September 2012 on the Arabian Peninsula. The virus causes severe and often lethal respiratory illness in humans with an unusually high fatality rate. The N-terminal domain (NTD) of receptor-binding S1 subunit of coronavirus spike (S) proteins can recognize a variety of host protein and mediates entry into human host cells. Blocking the entry by targeting the S1-NTD of the virus can facilitate the development of effective antiviral drug candidates against the pathogen. Therefore, the study has been designed to identify effective antiviral drug candidates against the MERS-CoV by targeting S1-NTD. Initially, a structure-based pharmacophore model (SBPM) to the active site (AS) cavity of the S1-NTD has been generated, followed by pharmacophore-based virtual screening of 11,295 natural compounds. Hits generated through the pharmacophore-based virtual screening have re-ranked by molecular docking and further evaluated through the ADMET properties. The compounds with the best ADME and toxicity properties have been retrieved, and a quantum mechanical (QM) based density-functional theory (DFT) has been performed to optimize the geometry of the selected compounds. Three optimized natural compounds, namely Taiwanhomoflavone B (Amb23604132), 2,3-Dihydrohinokiflavone (Amb23604659), and Sophoricoside (Amb1153724), have exhibited substantial docking energy >-9.00 kcal/mol, where analysis of frontier molecular orbital (FMO) theory found the low chemical reactivity correspondence to the bioactivity of the compounds. Molecular dynamics (MD) simulation confirmed the stability of the selected natural compound to the binding site of the protein. Additionally, molecular mechanics generalized born surface area (MM/GBSA) predicted the good value of binding free energies (ΔG bind) of the compounds to the desired protein. Convincingly, all the results support the potentiality of the selected compounds as natural antiviral candidates against the MERS-CoV S1-NTD.

Bioinformatics studies on a function of the SARS-CoV-2 spike glycoprotein as the binding of host sialic acid glycans.

SARS-CoV and SARS-CoV-2 do not appear to have functions of a hemagglutinin and neuraminidase. This is a mystery, because sugar binding activities appear essential to many other viruses including influenza and even most other coronaviruses in order to bind to and escape from the glycans (sugars, oligosaccharides or polysaccharides) characteristic of cell surfaces and saliva and mucin. The S1 N terminal Domains (S1-NTD) of the spike protein, largely responsible for the bulk of the characteristic knobs at the end of the spikes of SARS-CoV and SARS-CoV-2, are here predicted to be "hiding" sites for recognizing and binding glycans containing sialic acid. This may be important for infection and the ability of the virus to locate ACE2 as its known main host cell surface receptor, and if so it becomes a pharmaceutical target. It might even open up the possibility of an alternative receptor to ACE2. The prediction method developed, which uses amino acid residue sequence alone to predict domains or proteins that bind to sialic acids, is naïve, and will be advanced in future work. Nonetheless, it was surprising that such a very simple approach was so useful, and it can easily be reproduced in a very few lines of computer program to help make quick comparisons between SARS-CoV-2 sequences and to consider the effects of viral mutations.